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Abstract 

We establish from first principles a perturbative framework that allows us to compute reaction rates for processes 
taking place in nonequilibrium 0{N) linear-sigma systems in broken phase. The system of our concern is quasiuniform 
system near equilibrium or nonequilibrium quasistationary system. We employ the closed-time-path formalism and 
use the so-called gradient approximation. No further approximation is introduced. In the course of construction of 
the framework, we obtain the gap equation that determines the effective masses of tt and of a, and the generalized 
Boltzmann equation that describes the evolution of the number-density functions of tt and of a. 

ll.lO.Wx, ll.30.Qc, 11.30.Rd 

I. INTRODUCTION 

Lattice Quantum Chromodynamics (QCD) results indicate that chiral symmetry is spontaneously broken at T ~ 150 
MeV. Such temperatures may be reached in relativistic heavy-ion collisions. Thus, the time evolution of the chiral 
phase transition may be traced through observing predicted consequences, e.g., in the production rates of dilepton 
and of photon, etc. |]l|, or the formation of disoriented-chiral-condensate (DCC) domains (see, e.g., Q). 

In this paper, for an effective low-energy model for QCD, we take the 0{N) linear-sigma model. The case of = 4 
is of practical interest. The ultimate goal is to construct a framework for analyzing how the phase transition proceeds 
through DCC's. A numerous work was reported so far on the "static" properties of this model or its variants. (See, 
e.g., 1^,^. Earlier work is quoted therein.) Since Rajagopal and Wilczek ||] proposed that heavy-ion collisions can 
generate large DCC domains, analyses of dynamical evolution of the system from a symmetric phase to a broken 
phase have energized. The analysis in jsj has been refined in and has been extended j^] by incorporating more 
realistic descriptions of heavy- ion dynamics. The effects of quasiparticle excitations are studied in Refinements 
employing realistic initial conditions are made, e.g., in Since then, much work has been devoted to the analysis by 
systematically taking quantum and medium effects into account: For example, analysis on the basis of density-matrix 
formalism is made in large- A^ limit of the 0{N) linear-sigma model has been studied in closed-time-path 
(CTP) formalism of nonequilibrium dynamics is employed in [ p^ , Calderira-Leggett theory is applied in jist , a self- 
consistent- variational approach has been taken in , the time evolution of a particle distribution is studied in [[l5| , 
and relaxation rate for long- wavelength fluctuations are analyzed in (For other related works, see, e.g., |[17[.) 

Different assumptions and approximations are employed in these analyses. 

In this paper, as a first step toward the goal, we lay down from first principles a perturbative framework on the 
basis of a loop-expansion scheme. Only approximation we use is the so-called gradient approximation (see below). We 
use the standard framework of nonequilibrium statistical quantum-field theory that is formulated by employing the 
closed-time path, —oo — > +oo — > —oo, in a complex-time plane |l^-|20[| , which is referred to as the CTP formalism. 

Throughout this paper, we are interested in quasiuniform systems near equilibrium or nonequilibrium quasistation- 
ary systems. Such systems are characterized by two different spacetime scales; microscopic or quantum-field-theoretical 
and macroscopic or statistical. The first scale, the microscopic correlation scale, characterizes the reaction taking place 
in the system, while the second scale measures the relaxation of the system. For a weak coupling theory, in which we 
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are interested in this paper, the former scale is much smaller than the latter scale.fj In a derivative expansion with 
respect to macroscopic spacetime coordinates X, we use the gradient approximation throughout: 

F{X, ...) ~ F{Y, ...) + {X- YrdY.F{Y, ...) . (1.1) 

Let A{x,y) be a generic propagator. For the system of our concern, A{x,y), with x — y fixed, does not change 
appreciably in {x + y)/2. We refer the first term on the R.H.S. to as the leading part and the second term to as 
the gradient part. The self-energy part Ti{x^y) enjoys a similar property. Thus, we choose x — y as the microscopic 
coordinates while X = (x + y)/2 as the macroscopic coordinates. 

The plan of the paper is as follows. The 0{N) linear-sigma model is introduced in Sec. II and the forms of retarded 
and advanced bare propagators are given in Sec. III. In Sec. IV, a perturbative framework is constructed from first 
principles. The framework thus constructed allows us to compute reaction rates by using the reaction-rate formula 
pl[ . In Sec. V, a quasiparticle representation of the propagator is given. In Sec. VI, after constructing the self-energy- 
part resummed propagator, the gap equation and the generalized Boltzmann equation are derived. Section VII is 
devoted to conclusion and outlook. Concrete derivation of various formula used in the text is made in Appendices. 



II. 0{N) LINEAR-SIGMA MODEL 

The Lagrangian (density) of the 0{N) linear-sigma model reads 



\2 2 T2 



4! 



B ■ <PB 



(2.1) 



where (t)B — (^^j, 0^, 0^). When < 0, £ describes the system whose ground state is in a broken phase 
in the classical limit. Hb is an external field, which explicitly breaks 0{N) symmetry. Noticing the fact that a 
renormalization scheme for the symmetric phase (m^ > 0) works [ p2[ as it is for the broken phase (m^ < 0), we 
introduce renormalized quantities, 4>b = \fZ4'^ = Zmm^, and = Zx\, in terms of which C reads 



A 



2 22 



1 



r2\2 



(2.2) 



where H ^/ZHB■ 

The system in the broken phase is governed by the "sifted Lagrangian:'' 



C{(j,{x)-^{x)) = C{(j>{x) + 0{x)) - 



dC{0{x)) 
d0{x) 



To avoid too many notations, for C{(j){x);0{x)), we have used the same letter "£" as in Eqs. ( ^.1[ ) and (2.2). (p{x) 
is the (classical) condensate or order-parameter fields and 0(a;) is the quantum fields, which describes the fluctuation 
around ip{x). Straightforward manipulation yields 



C{(l){x);(p{x)) = a + Crc + 



(2.3) 



^ 2 



A 
12 



(^202 + 2(V-^TJ"^(V' 



Crc = \{Z- l){d$f - \{Z,^Z - l)mV 



-(ZaZ2 - 1) 



12 



3! 



4! 



(2.4) 



^ It should be noted, however, that, as the system approaches the critical point of the phase transition, the microscopic 
correlation scale diverges. Thus, the formalism developed in this paper applies to the systems away from the critical point. 
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In Eq. (2.3), stands for the terms that includes only c-number field (p, which is not necessary for the present 
purpose, but plays a role when spacetime evolution of the system is studied. We ignore throughout in the sequel. 
It is to be noted that C enjoys 0{N — 1) symmetry. 

For obtaining an efRcient or rather physically-sensible (perturbative) scheme, we introduce P,^-^ weakly x- 
depcndent masses, Mt^(x) and AIa-{x), where 'x' is macroscopic coordinates. A consistent perturbative scheme is 
obtained by assuming that, when compared to m^, x^i^) = M^(x) — = tt, a) are one-order higher in the loop 

expansion. How to determine M^(x) and M^{x) will be discussed in Sec. VIB. We then rewrite C{(j){x); (fix)) in the 
form, 

C{(j){x); ip{x)) = CoWx); ip{x)) + Lint + + C„^c ■ (2.5) 



Here Crc is as in Eq. (2.4) and 



jro{${x);0{x)) - i / d*y'${x)A~Hx,y)${y), (2.6) 



C.nt = -^C0-$)$' - - €^ + , (2.7) 

C,nc = ^'$ix) ixAxWAx) + xAx)PAx)] , (2.8) 



where 



A-\x, y)^- [{d' + Mlix)) P^x) + {d' + Mlix)) P.(x)] 6\x - y) , (2.9) 
PAx)=l-\^{x)){^{x)\, VAx)^\^{x)){^{x)\, (2.10) 

Ml{X)^Ml{X) + ^0^{X), Ml{X)^Ml{X) + ^0^{X). (2.11) 

In these equations, the boldface letters denote N x N matrices act on the real vector space, Ptt (Pa) is the projection 
operator onto the tt- (cr)-subspace, and I is an unit matrix. |i^(a;)) is a unit vector along (p{x), and {(p{x)\ is an adjoint 
of|<^(a;)). 

The construction of perturbation theory based on £, Eq. (|2.5|), starts with constructing the Fock space of the quanta 
described by (j), which is defined "on <^(a;)." As stated in Sec. I, we are concerned about the systems which are not 
far from equilibrium states or from stationary states. Then, the theory to be developed may be applied to the case 
where (p{x) changes slowly. In one word, 'x' of (p{x) is macroscopic coordinates. 



III. RETARDED AND ADVANCED PROPAGATORS 



For the purpose of later use, we construct retarded and advanced propagators. The retarded (advanced) propagator 
A/f (Ayi) is an inverse of A~^(x,?/), Eq. (2.9), under the retarded (advanced) boundary condition. Derivation of the 
form of A/j and A^ in the gradient approximation is straightforward and, here, we display the result: 



Afl(^)(x,2/) ~ P^{x)A^^I^^{x,y)PAy) 



A 



(1) 



-iP-{x-y) 



A 



A^\)^)(x,y) = 2z(|^(X)) {V{X)\) j ^e-^^<^-y^P,A''^\X;P)A^;^\X;P). 



(3.1) 
(3.2) 

(3.3) 



In Eq. (3^), summation over ^ runs over n and a, and, in Eqs. (3^) and (3^), X — {x + y)/2, dx=dx — dx, and 



A 



(?) 

R{A) 



iX;P) 



P^ - MliX)±ie{po)0+ 



That the gradient part A^^'^j (a;, y) appears is a reflection of the fact that the internal reference frame, which defines 
a mode and three- tt modes, at the spacetime point x is different from that at the point y. Let us see the meaning of 
\(p)P- d {0\ in A^^,''^^. From Eq. ( 2.10 ), we can easily see that 
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p,|^)p-a(^|p. = l^)p-5(^l, 
p^\^)p-d{0\Pc = o ((e,c)^(a,^)), 

where use has been made of tp ■ d(p — [if = (p/\ip\). \(p)P- d (<^| enjoys similar property. Then, \(p)P ■ d{Lp\ and 
\0)P- d (<^| "induce" transition between tt mode and a mode. 



IV. CONSTRUCTION OF PERTURBATIVE FRAMEWORK 



A. Preliminary 

The CTP formahsm is formulated p9| | by introducing an oriented closed-time path C (= Ci +C2) in a complex-time 
plane, that goes from —00 to +00 (CT) and then returns from +00 to —00 (C2). The real time formalism is achieved 
by doubling every degree of freedom, (jj — > {4>i,4>2) and Lp —> ((^1,(^2), where (/)i(xo,x) = (/)(a;o,x) with xq £ Ci and 
(j)2{xo,x.) = 0(a;o,x) with xq € C2, etc. A classical contour action is written in the form 



dxQ / dx£{(t){x);ip{x)) = 



+00 



dxo / dx C{x) . 



C = C{cj)i\0i) -£(02; (^2) . 



(4.1) 



C here is sometimes called a hat-Lagrangian. 

Taking £oi which corresponds to Cq in Eq. (2^), for a free hat-Lagrangian, we construct from first principles a 
perturbative framework. Throughout this paper, we do not deal with initial correlations (see, e.g., [p^). Following 
standard procedure, the four kind of propagators emerges: 



Arf(x,y) = 



T 



-zTr 



(0?(x)0f(y) 
4{y)cj,-i{x) p 



P 



^22^ (a;, y) 



-iTr 



T(0^(x)0f(2/)) p 



(4.2) 



where p is the density matrix, T is a time-ordering symbol, and T is an anti-time-ordering symbol. At the end of 
calculation we set $1 = (j)2 and (fii = 02 ll9[| . Let us introduce a matrix propagator A.{x,y), where the bold-face 
denotes, as above, the N x N matrix and the 'caret' denotes the 2x2 matrix: A"j^ is the (a, /3)-component of the 
N X N matrix and, at the same time, (i, j)-component of the 2x2 matrix A"''. The matrix self-energy part 

S(a:,2/) is defined similarly. 

We stress again that the argument 'x' of F^{x) (^ — tt, cr) in Eq. ( ^.9[ ) is macroscopic spacetime coordinates. Let 
A be a propagator or a self-energy part. Due to 0{N — 1) symmetry of C, Eq. ( |2.5| ), we may write A as 

A(x, y) P^(a;)i(«) (x, y)P^ (y) + T^^^(x)i|«) (x, y) , (4.3) 
where, for f = tt (ct), — a (tt) and 

T^,(a:) ^ \0{x))d.^ (^(x)| , T^J,(x) ^ ^(a;)) (^(a:)| (4.4) 
(cf. Sec. III). One can replace T^^(a;) in Eq. (4.3) with T^^(j/), since the arising difference is of higher order (cf. 



Eq. (LI)). Fourier transforming A'^^^ and with respect to x — y, we have 



k{x,y)^P^{x) J i^e-^-(--«)i(«)(X;F)Pay) + T^^(x) J i^e-^-(--^) (X; P) , 



(4.5) 



X-derivative and A^p {— Af {X; P)) contains explicit (first order) X^-derivative: 



where X = {x + y)/2. In general, vl^^^ in Eq. (L5) consists of two pieces, A'^^'^ = A'^^ + A''^\ where A'j^^ is free from 



A{x,y)^Pi{x) 



d^P 



(2^) 



4«)(X;P)P^(y)+P4(x) j ±f^e-^^<^-yn[^\X;P)Pdy) 



d-^P 
(2^ 



-'P-^'^-y) Af{x-P)^ 



(4.6) 
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The first term o n th e R.H.S. is the leading part of A while the second and third terms are the gradient parts (cf. 
above after Eq. (1.1 )). 



B. Propagator 



From the definition of A's, Eq. (4.2), with (pi — (j)2 (cf. above after Eq. (4.2)), we see that 



= 



(4.7) 



01 = 02 



holds. Then, out of four Ay = 1,2), three are independent, for which we choose 

Ah = All - Ai2 , Aa = All - A21 , Ac A12 + A21 . 

Setting (pi = 02 (= 4>)j have 

iAf (a;,y) = 0(a;o-z/o)Tr{[r(a;), cj^" iv)] p] ^ e{x^ - y^)[r{x). <P^{y)\ , 
zAf (x,y) = -e{y^ - xo)Tr { [0"(x), c^^iy)] p] = -0(yo - x^) [^(a;), ^^(y)] 
i^f{x,y) = -iTT[{r{x)<P^{y) + (Pf'{y)r{x)) p] • 



(4.8) 



Thus, Aj^ and are the retarded- and advanced-propagators, respectively, which have already been constructed in 
Sec. III. Ac is the correlation function. Expressing A in terms of them, we have 



A = - 



2 V Ai 



where 



A+ = 



-/^R + A, 



1 ±1 
±1 1 



^Aci+ 



(4.9) 



(4.10) 



From the definition of Ac and Eq. (4.2) with 0i = 02, it follows that 

(zAc(a;,y))*=zAc(a;,y), ^f{x,y) = /^P'^{y,x). (4.11) 
To leading order (of derivative expansion), it is known that Ac takes the form (cf., e.g., ||2^ ) 

Ac(x,2;) ^ V^{x) J i^e-^-(-^) (1 + 2/^^;^)) {AfiX;P) ~ AfiX;P)) P^{y) , 

where X = {x + y)/2 and, as will be seen below, is the real function that is related to the particle-number density. 
Then, to the gradient approximation, we may write 



Ac(x,y) = F^ix)Ai^\x,y)P^{y)+T^^^{x)Al^\x,y) , 



where, T'^^{x) is as in Eq. (4^) and 

A(«) = Ag) . (1 + 2U) - (1 + 2/^) . Af + Ag) 



(4.12) 



(4.13) 



where f stands for tt or cr, and A^^' is a gradient part. It is clear that, in Eq. (4.13), although f is a 'repeated index' 
on the R.H.S., summation should not be taken over ^. This type of equations appears frequently in the sequel. Here 
we have used the short-hand notation F ■ G, which is a function whose "(x, y)-component" is 

[F ■ G]{x, y)^ I d^z F{x, z)G{z, y) , (4.14) 
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with '1' in Eq. (4.13) the function whose (x, y)-component is S'^{x — y). For a given p, Ac is computed through 
Eq. (4.8). Eq. ( 4.12| ) with Eq. ( 4.13| ) is unders tood to be th e def ining equat ion o f and A^^-*. Physical meaning of 
is clarified later (see, Sec. VIC). Using Eq. (pl|) in Eq. ( p^) with Eq. (p^ , we obtain the relations: 

(Ai«) (x, y)) * = - {x,y), (a^?) (x, y)) * - - Ajf) {x,y) , 
A(«) (x, y) = aP (2/, x) , A(«) {x, y) = A^f ' (y, x) . 

Note that T^^A^^^ in Eq. (|1|) and A^^^^ in Eq. ( [4.13D are the gradient parts (cf. above after Eq. (HJ)). 
Applying f3(a2 + A^|)P^ (= f3P5(92 + 7\/{2)) to AcA+, we obtain 



where 



T3 (dl + Mlix)) P^{x)A,{x,y)A+ 



(27r)4 



-'i.P-{x~y) 



2t{f^,P^-Ml}AfiX;P) 



2i {P ■ dV^{X)) (l + 2/e(X;F)) {Af{X;P) - A^ {X : P) 



(P2 - A^|(X)) (a«)(X; P) + Tf^(X)A«)(X; P) 



+ : 



^2 



^+Ac(x,y)(aj, +7W|(y))P4(2/)f3 
d^p 



(27r) 



2z{/^(X;P), P2-M|(X)} a(«)(X;P) 

2^ (P . (X)) (1 + 2U{X- P)) (a(«) {X- P) - A^^ {X; P)) 
(P2 - X|(X)) A^f (X;P) - (p2 - A1|(X)) Tf^(X)A(5)(X;P) 



{f^{X;P),P'-Ml{X)} 



dU{X-P)d{P^-Ml{X)) dMX;P) djP^-MliX)) 



dX, 



2P-dxfdX;P) 



d.f^{X;P) dMliX) 



dXu 



dPt^ 



dX, 



(4.15) 



(4.16) 



Bare-N scheme 



The propagator matrix A is an inverse of — 7-3(8^ + M^)P^ (cf. Eq. ( |4~l| ) and Eq. ( plj ) with Eq. (2J)). Then, 
Eqs. (4.15) and ( |4.16D should vanish. From this condition, we obtain the following relations: 



{fi, P'-Ml}=0, 
(P2-X|(X))A(«(X;P) = 0, 
A(;)(X;P) = -A|:)(X;P) 



2iP„ 



(1 + 2UX; P)) ( a'^^\X; P) - A'^> {X; P) 



Ml{X)-Ml{X) 

- (1 + 2UX- P)) (a*,") (X; P) - A^")(X; P; 



(4.17) 
(4.18) 



(4.19) 



As will be shown later (cf. Eq. (6.8)), (^ — tt, cr) is related to the number density of ^: f^{X;TEp^\p) = 

—9{—t) + e{T)N^{X] Ep^\ e{T)p) (r = ±). Then, Eq. (4.17) is a "free Boltzmann equation." One can construct 
a perturbation theory in a similar manner as in p6| , where a complex-scalar field system with symmetric phase is 
treated. We call the perturbation theory thus constructed the bare-A^ scheme, since A'^ obeys the "free Boltzmann 
equation." This theory is equivalent [E6| to the one obtained in the physical- A^ scheme, to which we now turn. 
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Physical-N scheme 



We abandon Eq. (4.17), while we keep Eqs. (4.18) and (4.19). This means that in the present (physical- A^) scheme 
differs from in the bare-A'^ scheme. Specification of /c is postponed until Sec. VIC, where we require the number 
density to be as close as possible to the physical number density. Now, A is not an inverse of — T3(9^ + A^|)P^. 
It is straightforward to show in the gradient approximation that A is an inverse of —T^{d^ + Al|)P^ + L'(x,?/)y4_, 
where A- is as in Eq. (4.10) and 



L'(.,y)=zP,(X) j ^ 



Then the free action is 



1 



Ao = -- d^x(f>ix)T3 {d^+Ml{x)F^{x)) (b{x) + - / d4xd4y0(^)L'(x,y)A_0(y) 



(4.20) 



(4.21) 



Note that the Lagrangian density corresponding to the last term of Eq. (4.21) is nonlocal not only in 'space' but also 
in 'time.' Here it is worth mentioning the so-called |po|-prescription. With this prescription, at an intermediate stage, 
we have L'(a;, y), which is l ocal i n time. (For completeness, we briefly discuss the |po|-prescription in Appendix A.) 
Since the last term of Eq. ( 4.21 ) is absent in the original action, we should introduce the counter action to compensate 

it, 



Ac 



d^xd^y^{x)L'{x,y)A^Uy), 



(4.22) 



which yields a vertex —iL' A- (= iVc). From Eqs. 
propagator becomes 



and (4.22), Vc • A • Vc = follows, and the Vc-resummed 



^c—resum 



1 + E 



^V, • A 



A • V, • A 



A + SA 



Since Vc is a gradient part, we obtain, to the gradient approximation. 



dAc^iP^iX) 



(27r)4 



-iP-(x-y) 



{/,, P'-Ml}AfAfA^. 



(4.23) 



Note that Eq. ( 4.23 ) possesses pinch singularities in apo-planc, due to A^' {X; P)A''2' {X; P). Since SA is proportional 
to j4+, it contributes to Ac (cf. Eq. (4.9)). Then, including 6 A to Acj4+/2, we obtain for Ac (within the gradient 
approximation) , 

Ac(x,y) = -iPc(x) J ^e'^^<^-yh{po)il + 2f^(X;P))S{P'~MliX))P^{y) 



+ 



d^P 
(2^ 



(2^) 

-iP-(x-y) 



^iVdX){k, P^-Ml][Af{X-P)~Af{X-P) 



+P^(X)A(f (X;P) + {\^{X)) dx, mX)\) Al^HX;P) 



(4.24) 



where A]^' is as in Eq. ( p^ ). From Eq. ( [4. 18] ), we see that A^f^(X; P) cx S{P^ - MjiX)), so that the term with 
A'-£ in Eq. (mi) may be absorbed into the first term of the R.H.S. by modifying the definition of f^{X;P). Thus 
we shall drop the ter m w ith A^.^-* her eafter . 



(4.25) 



Substituting Eqs. (^, (^), and ( |4.24| ) into Eq. (|J), we obtain 

d'^P 



A{x,y) = A^''\x,y) + 
A(0)(x,y)=P,(x-) J 



(27r)4 

e-^P-{^-y)^(i){^X;P)I>^{y) 



A^P\X-P) + A^*\X-P) 



A(rt(X; P) = -^Pe(X) {/,, P2 - Ml) [Af - Af)"A+ , 
A^'\X-P) = 2i\^)P-lj 



(4.26) 
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where A^^^ = A'^'iX; P), = /^(X; P), etc., and 



A(«)(X;P) 



n = 



A 



A(A<,^^-A^^) 



(l + /e)(A(«)-A 



V A^-^A^,-) 
1 

Ml -Ml 



fa 



(0 



k (a«^ - A, 



(4.27) 



As has been observed in Sec. Ill, \'^)P ■ d{if>\ and \'f>)P- d (<^| in the gradient part A^*) "induce" transition between 
TT mode and a mode. Thus, A*^*' contributes, e.g., to the one-loop two-point function with one-7r and one-cr legs. 



C. Feynman rules 



Two fundamental elements of Feynman rule s are the propagator and the vertices, which take th e (2 x 2) matrix form . 
The propagator- matrix is given by Eq. ( 4.25 ). The vertex factors may be reads off from Eq. (41) with Eqs. (^.5|), 
(2.4), ( p7^ ) and ( |2.8| ). As a matter of fact, from Eq. (4.1), it is obvious that there is no vertex that mixes ipis with 
(/)2's. The vertex factors for the fields 0i's are the same as in vacuum theory. Each vertex factor for (/)2's is of opposite 
sign to the corresponding vertex factor for ^I's. Thus, in matrix notation as for the propagator, every vertex-matrix 
V is diagonal, V = diag.(v, — v), with v the same as in vacuum theory. All other elements of Feynman rules, e.g., 
integration over every loop momentum, are the same as in vacuum theory. 

Having thus constructed Feynman rules, one can compute reaction rates for processes taking place in the system, 
by using the reaction-rate formula pl| . 



V. QUASIPARTICLE REPRESENTATION OF THE PROPAGATOR 

Here we obtain a quasiparticle representation [ p7[ of A, which tremendously simplifies the practical computation. 
Straightforward but lengthy calculation shows that A may be written in the form (cf. the definition (4.14)): 



B, 



^R 









(5.1) 



where 



Bfl 



B,f + \^{X)) dx. (^(X)K(X) 
Bf + \cp{X)) dx. mX)\a>^{X) 



1 

1 

1 1 





with 



(0) 



B 



1 h 
1 i + 

2i 



B 



(0) 



1 + /? h 
1 1 



d'^P 



MUX) -MUX) J (2^)4 



e-^P-(-v)p^(f^{X;P)-fUX;P)) 



It can readily be shown that 



(5.2) 



Bl • tsBr ~ B_R • tsBl ~ Ta . 



(5.3) 



It is obvious from the argument in Appendix A that, at an intermediate stage before taking the |po|-prescription, 
Bi and Bij are local in time and satisfy the relation (5.3). Then, B^ and Bij at that stage are the generalized 
Bogoliubov-matrices 12^. Eq. (5.1) tells us that, through Bogoliubov-transforming (fi {i — 1,2), one can introduce 



8 



the stable quasiparticle modes, 0J {i = 1,2), whose matrix-propagator is diag.(A/j, — A^i) in Eq. (5.1) (cf. for 
more details). 

It is worth mentioning that and may be written as 



where 



Bi(x,y)~P^(2;) 



S'^{x - y) 



S {x~y) S*{x~y) + F^{x,y) 



P?(y): 



Br{x. 



y)-fi[x)y ,54(2; ~y) S\x-y) J^^^y'' 



dxdx fd^.y) + h{x,y) dydy 
Ml{X)-Ml{X) 



with TT = (T and a — t:. 



VI. THE GAP EQUATION AND THE GENERALIZED BOLTZMANN EQUATION 
A. Self-energy-part resummed propagator 



We write the bare propagator A in Eq. (|^) as A = A(°) + A^^) (A^^) = A^p) + A^*)). A^") is the leading 
part and the gradient part A*^^^ represents variation in the macroscopic spacetime coordinates X^, through first-order 
derivative dx^ ■ Interactions among the fields give rise to reactions taking place in a system, which, in turn, causes a 
nontrivial change in the number density of quasiparticles. Thus, the self-energy part S ties with A'"'^-'. More precisely, 
is of the same order of magnitude as A^^^. Hence, in computing S in the approximation under consideration, it 
is sufficient to keep the leading part (i.e., the part with no X^-derivative) : 



tix,y) Pdx)m{x,y)Pdy) P^X) J 



(2^ 



-y)t(i){X;P) . 



(6.1) 



We are adopting the loop-expansion. Then, to the gradient approximation, the relevant self-energy diagrams are 
the one-loop diagrams, together with relevant counter diagrams. There are two one-loop diagrams, the one is the 



tadpole diagram that includes one vertex coming from c\'^f (cf. the definition (2.7)), and the one is the diagram 



^ ("3") 

with two v erti ces co ming from C^J^. The counter diagrams are the diagrams that include £rc + ^rc (cf. Eq. (2.5) 



with Eqs. ( |2.4| ) and (|2.8[ ).) The contribution from the tadpole diagram includes no X^-derivative. Yet higher-order 
contributions to S come in when one proceeds beyond the gradient approximation. 
A S-resummed propagator G obeys the Schwinger-Dyson equation: 



G = A + A S G. 

In Appendix B, we solve this within the gradient approximation. The result is 

d^P 



(6.2) 



G(x,y)~G(0)(a;,y) 



(2^)4 



G^P\X;P) + G^'\X;P) 



P^x 



d^P 



27r)4 



<^-y^&^\X;P)Pdy), 



GY>{X;P) + G\P>{X;P). 



-zP^X) 



G 



A, 



G(o)(x,y) 
&P\X;P) 
g['\X;P) 

G^i\x-P) 

G(*)(X;P) = 2z|^)P. d {v\ ( ^w^T-f nMMa) 1 +2* W)P ■ d{^\uG ~ ^P- d {^Wg 

V ^ '^A ^^A '^A 



}( 



P'-Ml-Re^f}] GfGfA^^ 



(7r)^((T) 







(6.3) 
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where 



G 



(0 



G 



(0 



zr^^ = (i + /,)siV-/eS^V, 



G 



(?) 



/? (Gk^^ 



G 

(?) 



G 



-G 



G' 



In the above equations, = /^(X; P), E's = P)'s, A^^ ^ MjiX), and 



G 



(?) 



G 



(?) 



v(«) - v(«) 
^fl;(A) ~ ^11 



-(?) 



12(21) 



In Appendix C (cf. Eq. (|C5D), we show that T.f{X; P) = (^^^^^(X; P)j . The expression for Tr^ is given in Appendix 
D, and E^l^ and S^^^^ in Eq. (jej) on the mass-shell are computed in Appendix E, which plays a role in Sec. VIC. 



P2-M|(X)-4%(X;P) 



(6.4) 
(6.5) 



(6.6) 



.(?) 



B. Gap equation 

We have introduced two "mass functions" Af|(X) (^ = t^,o'), which is a generalization of [Q, in which only one 
mass parameter is introduced for computing an effective potential for the equilibrium system. In this subsection, we 
determine so far arbitrary mass functions. Various methods are available to this end (see, e.g., Q ). Among those we 
employ the on-shell renormalization scheme. 



ReY.f{X;pt = Ml{X), p = 0) = 



s.p. 



(6.7) 



(^ — TT, (t). For a subtraction point, 's.p.,' instead of (pq — A^?(A"),p 



0) adopted here, we can choose [p^ — 
p1 + A^|(Ar) with Ps arbitrary. Eq. ( |6.7[ ) is the gap equation, by which M|(A') 
(f = TT, cr) are determined self consistently. The explicit form of the gap equation (|6.7D to leading one-loop order is 



Ps), where E'^f 



displayed in Appendix D. M^{X) thus determined depends on Lp^{X) and /t. With M^{X) in hand, we can judge 
if A4|(A'), Eq. ( 2.11 ), is positive or negative. If A^|(A') < 0, the log-wave-length modes, < \M({X)'^\, have 
imaginary frequencies and thus unstable, which causes perturbative instability of the system psf . Unfortunately, no 
consistent scheme for treating such a case is available. Then, when A4'^{X) < happens, we abandon the condition 
(3.7) and set Ai'^{X) = 0. Perturbation theory in this case is less efficient. 



C. Boltzmann equation 

In order to find a physical meaning of , we recall a momentum density of the system: 



Vix) ^ -Tr 



d(l){x) 
dxo 



V0(y) } p 



J y=x 



i__d_ 

2 dxQ 



y = X 



We first analyze the contribution from G^c^ {— G^^^ + G)22)i Eq. ( |6.3[ ). When the interaction is switched off, G, 



i(0) _L n(o) 



reduces to A^"' (— A}-^^ + AjgOj Eq. (|4.26|). Computation of the contribution from Ac"' yields 



,(0) , a(o)^ 



.(0) 



(0) 



nx) 



C(, f d^p 

A<") ~ T y (2^ 



^p{[/^(x;4«),p) + 1/2] + [/^(x-;-4«),p) + 1/2]} 
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= 1, = 



where c,r = A'^ — 1, c, 
free (A = 0) energy density yields 



and p = p/|p|. Similarly, the computation of the contribution to the 



Subtracting the contribution from the vacuum, we see that 
of ^ through 



TT, a) is related to the number density N(^{x; eI^\ p) 



N^{x, Ep,p) = f^ix; p) = -l-f^ix; -E^^-p) . 



(6.8) 



It is to be noted that the argument 'x' here is mac rosc opic coordinates. For the interacting system, the corresponding 
relation is obtained using Gc (= Gn + G22), Eq. ( |6.3D , and the total energy density V^ix) in place of Vj^^^ix) above. 



The contribution to V^{x) from the difference G 



(0) 



.(0) 



yields a correction to the relation mM). This is also the 



Since (G^^^),cxGi«^G;4 , 



case for the contribution of {G-^ )c- (G^*^)c yields vanishing contribution to the gradient approximation. 

We now turn to analyzing the remaining contribution that comes from(G^^^)cinEq. (fj 
in the narrow- width approximation, IttiTS^ [X] P) — )■ — e(po)0^, pinch singularity is developed. Then, the contribution 

of {G'f^)c to T"^{x) diverges in this approximation. In practice, JmE^^ (cx A^) is a small quantity, so that the 
contribution, although not divergent, is large. This invalidates the perturbative scheme and a sort of "renormalization" 
is necessary for the number density J 



This observation leads us to introduce a condition {G'f'')c = or 



(6.9) 



This serves as determining equation for so far arbitrary /j. Then G in Eq. ( |6.3| ) becomes G — G*^"^ + G^^' + G'^*^ 
which is free from pinch singularity in the narrow- width approximation. It is obvious that, in the present scheme, 
above-mentioned large contributions do not appear. 

In order to disclose the physical meaning of Eq. (3.£), we first define on the mass-shell, po = ±(jj^(X; ±p) (= ±^^5^): 



P'^-M'I{X) - ReT.^^\X;P 



= 0. 



(6.10) 



From Eq. (3/?), we see that, when M'j{X) > 0, uj'£^{X; 0) = M^{X). We also introduce a wave-function renormal- 
ization factor, 



V^i- 



1 



dReT. 



(0 



Po = i^iiX- p) 



2uj^{X;p) dpo 

It is now straightforward to show that Eq. ( |6.9[ ) becomes, on the mass-shell, 
dNf _ _ duj^{X;p) dN^ 



dXn 



dXn ^P^' 



Po = i^aiX; p) 



dN^{X:uj^{X;p),p) duj^{X;p) dN^ duj^{X;p) dN^ 



dXn 



Po = ^diX; p) 



<P) 



Po = ujs(X;p) 2a;^(^;p) 



i9p 



ax 



ax 



dp 



(«) 
12 



Po = uJi^(X; p) 



(6.11) 
(6.12) 



Here is as in Eq. ( |6.8| ) with Ep^' Lu^{X;p) and = duj^{X;p)/dp is the velocity of the quasiparticle mode 
with momentum p. F^^^ in Eq. (6.12) is the net production rate of the quasiparticle of momentum p. In fact, F^^' 



is the difference between the production rate and the decay rate, so that F^ is the net production rate. In the case 
of an equilibrium system, F^^-* — (detailed balance formula). iV^ — N^(X;u;^{X]p),p) here is essentially (the main 



part of) the relativistic Wigncr function, and Eq. (6.11) is the generalized relativistic Boltzmann equation (cf. |29|). 
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Let us suppose the case >> i?eS^M The sohition to Eq. (|6.10|) is po ~ ±Ep^\ To one-loop 



order under consideration, only diagram that contributes to s|^2|2i) ^'i- ( 6-12 ) is the one that includes two vertices 
coming from Thus, ^^2(21) contains two Ai2(2i)'s, each of which contains on-sheU (5-function, 5{Q^ — A4'^{X)) 

(cf. Eq. ( 1.27] )). One can easily see then that I]^2}2i) I ^ ^(j, vanishes unless Alcr > 27U7r- Tif^f^^-^^lX; P) at po = Ep^'' 
— TT, (t) is computed in Appendix E. Let us turn to the case where M'^{X) < O^S^"*! ^ (j,^ with ^ = tt or 
^ = (T or ^ = TT and — a. For a hard momentum, p > O^E^^j _^ ^^-|^, the same statement as above holds. 

For computing '^'"12(21) i-^'^ ^) ^i*'^ ^^^^ boL P < _ j. («) ) ' above two Ai2(2i)'s contained in S^^Li/s 

should be replaced by Gi2(2i)'s [|o|. 

What we have shown here is that the requirement of the absence of Gj^'* oc G^'*G^'' from G leads to the Boltzmann 
equation for the quasiparticle-distribution functions. This means that the quasiparticles thus defined are the well- 
defined modes in the medium, in the sense that no large contribution appears in perturbation theory. Conversely, if 
we start with defining the quasiparticles such that their distribution functions subject to the Boltzmann equation, 
then, on the basis of them, well-defined perturbation theory may be constructed. 

Comparison of our derivation of the generalized Boltzmann equation (GBE) with those in related works is made 
in pGf . It is worth recapitulating here the comparison with the derivation in nonequilibrium thermo field d ynam ics 
(NETFD) which is a variant of nonequilibrium quantum field theory. We have imposed the condition (6^) for 
determining or the number density N^, so that pinch singularities (in narrow- width approximation) disappear. On 
the other hand, in NETFD, which employs the (space)time representation, the GBE is derived by imposing "an on- 
shell renormalization condition" for the propagator. Since the pinch singularity is a singularity in momentum space, it 
is not immediately obvious how to translate this condition into the (space)time representation, as adopted in NETFD. 
Nevertheless, closer inspection of the structure of both formalisms tells us that our condition is in accord with the 
on-shell renormalization condition in NETFD. Incidentally, reconciliation of the NETFD with the |po I -prescription 
(cf. Appendix A), a notion in momentum space, remains as an open problem. 



VII. CONCLUSION AND OUTLOOK 



In this paper we have constructed from first principles a perturbative framework for computing reaction rates of 
the processes taking place in the 0{N) linear-sigma system in a broken phase. Only approximation we have employed 
is the so-called gradient approximation, so that the framework applies to the quasiuniform systems near equilibrium 
or the nonequilibrium quasistationary systems. 

The reactions taking place in the system causes a spacetime evolution of the system — development of phase 
transition. This is the next subject following to the present analysis. At the final stage of such an analysis, one should 
check whether or not the rate of the phase-change of the system is too large, so that the gradient approximation 
adopted in this paper is violated. 
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APPENDIX A: IPol-PRESCRIPTION 



One starts with f^{x,y) that is local in time, f^{x,y) = S(xq — yo) X]t=± ff^^H^i ^^o)- Here ± denotes posi- 



tive/negative frequency part and the time coordinate 'xq' is of macroscopic. Then, in place of Eq. (4.20), we have, 
with obvious notation, 

L'{x, y) 0, iP^iX) J ^ e'P<-~y^ [m' {xo - yo)dx„9l^\X; p) 
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+S{xo - yo) 2p • 



dMliX) 



(Al) 



with Xq = xq. The form ( JAID is local in time. It is well known through the analyses of equilibrium case that the 
following replacement (|po [-prescription) should be made: 

with f^{X;P) as in ( 4.20| ) in the text. Adopting this prescription, we attain Eq. ( 4.2C| ). (See for more details.) 



APPENDIX B: DERIVATION OF EQ. (6.3) 



Here we solve the Schwinger-Dyson equation (6.2). Multiplying B^"^ ^, Eq. (5.2), from left and B)^^ ^ from right, 
we obtain 



,(o)-i 



where 



G = A + A S G, 



t^B^n^-t- Bf 



Straightforward manipulation using Eqs. (4.25) and (4.26) yields 



(Bl) 



-A. 



(B2) 



,Mx. y) ^ n^{X))dx. (^(A)| j ^ e-*^-(--^) (MX; P) - /.(A; P)) A^^\X; P)A^j\X; P) 

+2^|^(A)) (^(A)| j ^e-^^<^-y^P^iUX;P)-MX;P))A'^^\X;P)A^X\X;P) 



d'^P 



{k,P^-Ml]Af{X-P)Af{X-P) 



Using the identity ( pl| ) in Appendix C, we obtain 
-S^ 

Sfl — Sii + S12 , — Sii + S21 = — S22 — S12 , 

Soff = S12 • (1 + • /e • Pc) - Pe • • Pc • S21 + Sn • P^ • ./^ • p« - p^ • • P^ • Sn . 



(B3) 
(B4) 



Substituting the leading-order expression (6.1) for S's and using Eq. (CI) in Appendix C, we obtain, to the gradient 
approximation, 



Soff ~Pc(a;) 



-iP-{x-y) 



(2^)4 

d^P 



(H-/5(A;P))E(f (X;P) - /^(X; P)s(f (A; P) P^(y) 



2^^^"^/ (2.)4 



-){/,, I](«'-4f}p,(y), 



where use has been made of P^(aPj)P^ = (no summation over ^). To leading order under consideration, S^|-'(A; P) 
and S^f (A;P) are pure imaginary. Then, from Eq. (01, i?eE^f (A;F) = -i?ei;^^i^(A; P) and, from Eq. (p|) and 



Eq. (g), we obtain /mS^f (A;P) = /toI]^|^(A;P) = i(E^|^(A; P) + E^f (A; P))/2. Using these relations, we finally 
obtain 



Soff(:z:,y)c.P4(x) 1 ^^e'^^^'^-y) ^^^{X-P)V^{y) , 
I](«f'(A;P) ^ zf('')(A;P) + Pe4«)} . 



(B5) 
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Note that the first te rm o f the R.H.S. of Eq. (B5), being proportional to A^, is proportional to the net-production rate 
(cf. above after Eq. which causes the change in the number density. The second term is proportional to A and 

includes derivatives with respect to X^, and is of higher order. Then one can drop the second term. Nevertheless, we 
shall keep it in the following. 



Substituting Eqs. (B2) and ( p33| ) into Eq. (Bl) and retaining up to the terms that are linear in Soff, we obtain 

oo 

G = A + ^ A [-s • A " 







-S4 / — 



Gr 








Soff 




Gr 






Ga 



(B6) 



where higher-order terms have been dropped and 



'R{A) 



^R{A) 



+ A 



R(A) E {■'^R{A) ■ ^R{A)) 



Keeping terms linear in Aq// (Eq. (B2)), we can solve Eq. (B6) 
(G(a;,y))2i ==0, 

+2i\ip{X)) l3x^ {^{X)\ 



(2^ 



(G(a;,y))22 
(G(a;,y))i2 



R^ A 



2^\^{X))^x. mX)\ j ^ e-^-(--^) P^ (/.(X; P) ~ MX; P)) G^^^ {X; P)G^J^ (X; P) 



+2i\ip{X)) dx^ {^iX)\ 
d^P 



d^P 



-iP- 



P) -U{X;P)) g'^^\X; P)G^a\X; P) 



-zPe(X) 



(2^) 



rip) _ 



where G^^I^^^{X; P) is as in Eq. in the text. Computing G = B^''' • G • B^''' in the gradient approximation, we 
obtain Eqs. (|6.3D - (|6.5D in the text. 



Gf{X;P)Gf{X;P), 



,(0) 



APPENDIX C: PROPERTIES OF THE SELF-ENERGY PART 



It is obvious that the relation (4.7) holds for the full propagator G and also for the self-energy-part inserted 
propagator A • S • A. From the latter, one can readily obtain the relation: 



0. 



(CI) 



The expression for full Gr and G^ with (pi — 02 (= 4>) are given by Eq. (4.8) with Heisenberg fields for (j)'s. From 
Eq. (4.8), we obtain 



(a^J^) {x, y)) * = A^^) {x, v) , A^"(y, x) = A^ (x, y) 



(C2) 



Eq. (C2) is also valid for self-energy-part-inserted propagator, A • S • A 

\^R(A) ■ ^R(A) ■ ^R{A) ) - ^R{y 
. v/^'a' . \a'a\ /„ ^\ _ [ao 



^R{A) 


^R{A) ' 


^R 
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which yields 



a/3/3' xa'i 



Applying d + A^^P^ from both sides of Eq. (C3), we obtain 



(C3) 



(C4) 



As discussed in Sec. VIA, it is sufficient to compute the leading part of , Eq. (6.1). Then, from Eq. (C4), we obtain 

E^«'(X;P)^(4«'(X;P))* . (C5) 



APPENDIX D: ONE-LOOP FORMULA FOR Sh AND THE GAP EQUATION 

A. Form for Sr{X;P) 



He re w e display the concrete form of Hj^,{X:P) to one-loop order. For the relevant diagrams, see above after 
Eq. (3.1). For a UV-renormalization scheme, as in [Q, we use the MS scheme. Computation is a straightforward 
generalization of [Q and the final result reads 



Here 



Sfl(X; P) = P^(X)4«)(X; P) - xdX) , 



4-)(X;P) 



E^^^X; P)='- [{N - 1)1^ + 3I„] + '—^''{X) [{N - l)J^AX; P) + 9J,,(X; P)] 



J [{N + 1)1^+1,] + ^^\X)J^„{X- P) , 
6 9 



A 



A 
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J^TTa{X; P) 



MjiX) MjiX) 
In 

(47r)2 e^^l 

1 



d^q N^iX;EfiX),ci) 



(2^)3 



Ei^\x) 



(^ = TT, cr) , 



(4^) 



MUX) _ P^-MUX)+Ml{X) MUX) 



2P2 



MUX) 



S {S-P')'-{MUX)-MUX))' 
2P2 



iS + PU'-iMUX)-MUX))' 



(27r)3 



1 



E, 



NUX;E!j^\q) 



P^+Ml-Ml~ 2poE, 
NUX;E!j"\q) 



P^ + Ml-Ml + 2po4"^ - 2p • q + leipo + E\''>)Q+ 



2p-q + ie(po-4"^)0+ 



+ {(7 ^tt) 



J^^{X-P) = J^UX;P) 



M± = Mi ^ Mi 



(C = 7r,Cr) ; 



where Ejf\x) = Jq'^+M^iX) 



S 



{MUX) - MUX)f - P2 {MUX) + MUX)f - P2 



(Dl) 



(D2) 
(D3) 



and fid is an arbitrary para met er that appears in the dimensional-regularization scheme, as adopted here. The first 
term on the R.H.S. of Eq. (D2) is valid in the region P^ < 0. Its expressions in other regions of P^ are obtained 



through analytic continuation with Ml{X) M^{X) — ie{po)0^ 
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B. Gap equation 



The gap equation (S.7) with Eq. (pif) yields 



M?{X)-m'^ 



with 



Here 



and 



SMiiX) - TT, a) 



\ MUX) ^^MliX) 



{D4) 



+nii\X;po^M^iX),0) 



6 



[{N - 1)X„ + 



18 



{An)^\2MUX) MUX)' 
MUX) 



hi 



o2,,2 



M5 



(D5) 



(47r)2 



(iV - 1) In 



- ICiZ} ] + sVSn + 9 In 



MUX)] 



{(iV - l)nii\X;po = X.(X),0) + 9Hl'3j(X;po = A^.(X),0)} 



(D6) 



. . (A.^(X) - ,MliX)) VA^^(^)(^^W-4A^^(^) ) 



In 



- ^UX){MUX) - (j^)) 
" MUX) + ^MUX){MUX)^iMUX)) 



hi 



-2MUX) - ^MUX)iMUX) - ^MUxJ) 
MUX) ^2MUX) + VMUX)iMUX) - 4MJ(X)) 



(amUx) ~ MUX)) 

e{MUX)) arctan 



v/.M^(X)(4>(^(X)-Al^(X)) 
MUX) 

^MUX){iMUX)-MUX)) 



IC'-° 



-e{MUX) - ^MUX)) arctan 

-7:e{2MUX)~MUX))] , 
' {MUX) - 4MUX)) 



MUX) 

^MUX){4MUX)-MUX)) 
MUX)~2MUX) 



^MUX)(MUX) ^MUX)) MUX) ~ ^MUX)(MUX) ^MUX)) 



MUX) 

+e {4mUx)-mUx)) 



MUX) + VMUX) iMUX) 4MUX)) 

^MUX){iMUX)-MUX)) 
MUX) 



2 arctan 



^MUX)i'iMUX)-MUX)) 



MUX) 



Hf;(X;po,0) = 



(2^ 



MUX)NUX;Ei^\x)A) 



E, 



4^\x){mux)-4pUe!,^\x))^} 

{2MUX) - MUX)) NUX; E^^\x), q) 

{{2MUX) - MUX))' - 4pg(4'^)(X))2} 



Hjf (X;po,0) = 2 



d^q Ml{X)N^{X;Ei^\x),q) 
(2^)' Ei^\x) {mux) 4p§(4«)(X))2} 
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APPENDIX E: COMPUTATION OF S12 AND S21 



Here, we compute the leading part of the one-loop contribution to Si2(2i): 

Si2(2i) [X; P) = P^{X)J:[ll^^ {X;P) + P,(X)s(^j2^) {X- P) , 

J ^12(21) (.^'yi^i2(2i)(.^'-^^yi 



^12(21) (^'-P) 



^12(21) (^'-^) 



-^^^ / [(A^-i)AS2r)(^;Q)AS2,(^;^-Q) 

+9A[l\^,^iX;Q)A[l\^^^iX;P^Q)_ . 



We compute S^2'(2i) mass-shell po = Ep^\ The contribution to ^^2(21) from the term that accompanies two 



(2±^ ^"j- ±^[^^± f 

^^^^^^ ^ ^ only when 2A^^ < Ada- Computation is straightforward but 

lengthy. We only display the final forms: 



^^u\x;P) = 



727 



pq=Zi 



{X; P) = ^ / " N^iX; ^, q) [l + N^iX; ^ + E'^^ , qTp) 

727r p J^^^ L 



4l\X;P) 
4l\X;P) 



IUtt 

^44^^ m; — ^ 

144^ P 7^,3 



p ■ q = Z3 



2 + iV4X;4-)-C,p-q) 



p ■ q = Z3 



(El) 



with 



Ul(il) = 2Ml 



C«2(i2) — iul(ll) — Ej 



MgE'-^'^ ±-p^ MI-AMI 



2Mo 



Z2 



Z3 = 



2E^^\-Ml 

2p^e - Ml ' 

2E^^k^_Mlj^2Ml 

2p^e -Ml 

2e\^\-MI 

2p^e ~ Ml ■ 
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